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Abstract The beam energy is measured in et e~ collision by using Compton Backscattering. The 
uncertainty of this measurement process is studied by virtue of analytical formulas, and the special 
effects of variant energy spread and energy drift on the systematic uncertainty estimation are 
also studied with Monte Carlo sampling technique. These quantitative conclusions are especially 
important for understanding the uncertainty of beam energy measurement system. 
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1 Introduction 


The upgraded Beijing electron-positron collider (BEPCII) is a 7-charm factory with a center 
mass of energy ranging from 2.0 to 4.6 GeV and a design peak luminosity of 10?? cm~? s^! 
at beam energy of 1.89 GeV [1, 2]. The upgraded Beijing spectrometer detector (BESIII) with 
high efficiency and good resolution for both charged and neutral particles was constructed and 
started data taking in 2008 [3]. BESIII research region covers charm physics, charmonium 
physics, spectroscopy of light hadrons and 7-lepton physics [4]. 

After vast amounts of data are acquired and analyzed, the statistical uncertainties in analyses 
of physics become smaller and smaller, while the systematic uncertainties play more and more 
prominent roles [5, 6, 7], one of which is the uncertainty due to the measurement of beam 
energy. To reduce such an uncertainty, start from year 2007, a high accuracy beam energy 
measurement system (BEMS) was designed, constructed, and put into operation at the end 
of year 2010 [8, 9, 10, 11], which is of great importance for many physics analyses at BESIII, 
such as 7 mass measurement, charmonium resonance scans, and determination of the branching 
ratio with the uncertainty at the level of 1% to 296. The measurement procedure of BEMS 
can be recapitulated as follows [12]: firstly, the laser source provides the laser beam and the 
optics system focuses the laser beam and guides it to make head-on collisions with the electron 


(or positron) beam in the vacuum pipe, after that the backscattering high energy photon will 
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be detected by the High Purity Germanium (HPGe) detector, which is the key instrument of 
BEMS. The accuracy of beam energy depends merely on the detection result of backscattering 
photon. 

The essence of working principle of BEMS is Compton backscattering process (CBS). In 
order to understand the main feature of BEMS, the uncertainty of this measurement process 
is addressed by virtue of analytical formulas, where some experimentally meticulous details 
are neglected. These acquired quantitative results are of greatest consequence for the quali- 
tatively understanding the actual uncertainty of BEMS. Moreover, an experimentally special 
phenomenon is studied by simulation approach, which reveals a possible source of systematic 
uncertainty of BEMS. 


2 Energy formulas 


Here considered is a special and crucial case of CBS, that is the electron makes a head-on collision 
with the photon, whose geometry is sketched in Fig. 1. The energies of electron and photon 
before (denoted by subscript 1) and after (denoted by subscript 2) the collision are denoted as 


€1,2 (for electron) and w 2 (for photon), respectively. 
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Figure 1: Geometry of electron (e) and photon (y) before (a) and after (b) the head-on collision. 
The stright line denotes e while the wavy line y. 


In the light of the special theory of relativity, the energy and momentum can be expressed 
as 
w 


for photon and 


MeC EUV 
ps $ Pe = -3 > (2) 


for electron. In above equations, h is Plant constant and Me the electron mass. With Eqs. (1) 


and (2), it is readily to obtain the kinematic for electron and photon collision system. According 


to the law of energy and momentum conservation, 


Wy + E1 =WQ+€2, (3) 
and 
Wy £1U] W2 £2U2 
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Based on the Eqs. (3) and (4), it can be obtained with simple algebra 
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For BESIII, the optimal energy point is at 1.89 GeV, the velocity of electron with such high 
energy is very closely to that of light (c), in another word, v1/c ~ 1 with the negligible error. 
With such an approximation, Eq. (5) is recasted as 
€ 
w = ey (6) 
e 
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In BEMS, w; is provdied by laser and wz is measured by HPGe detector, and £1 is the beam 
energy that is to be determined with high accuracy. From Eq. (5), it is worked out 
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3 Uncertainty formulas 


The start point of uncertainty analysis of this section is the two formulas obtained in the previous 
section, viz. Eqs. (7) and (8). For brevity, in this section we adopt the nature unity where c = 1 


(the subscript e is also suppressed for electron mass) and begin with the comparatively simple 


case, that is Eq. (8), by virtue of which it is immediately obtained 
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On the strength of Eq. (8), Eq. (9) can be rewritten in the more concise forms as follows 
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'The compact expression for uncertainty evaluation is as follows: 
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with factors f+ defined as 
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Now we turn to Eq. (7). The similar process as that for Eq. (8) could lead to the fairly 
cumbersome derivative expressions of £1 with respect to w2, w1, or m, which have been degraded 
into appendix A. Herein we present another recipe. Comparing Eqs. (7) and (8), it is clear that 
the Eq.(8) is just the first term of Eq.(7), therefore it is nature to find an uncertainty expression 
for Eq.(7) which could incorporates the result acquired based on Eq.(8). To this end, we return 
to Eq. (5). If a function of £1, g(£1), is introduced, Eq. (5) becomes 
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with 
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Noticing the similarity between Eqs. (8) and (15), it is immediately obtained 
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Next, from Eq. (16), it also easy to get 
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4 Statistical and systematic uncertainties 


(14) 


(17) 


Analytical formulas acquired in the previous section are the foundation for the uncertainty anal- 


ysis relevant to CBS. In principle, there are two types of uncertainty: statistical and systematic. 


In a nutshell, the statistical are those types of uncertainties that have a random spread, and 


their uncertainties decrease with augment of data; the systematic include everything else. In 


practice, it is not always easy to distinguish two types of uncertainty, and sometimes it is rather 


difficult to identify the feature of systematic uncertainty. For example, as to the each term 


in error formulas such as Eqs. (11) and (18), it could include both statistical and systematic 


uncertainties. From a pragmatistic point of view, we lay stress on the relative magnitude of 


each term instead of focusing on its feature, and try to figure out the leading contribution for 


the uncertainty evaluation. 


Table 1: Some parameters related to BCS. 


Parameter Central Value Relative error Reference 


(value scope) 


Me 0.51099828 MeV 2.153x10-8 [13] 
ul 0.114426901 eV — 8.739x10-? [14, 10] 
w2 (2-7) MeV 5 x10758 [12, 10] 


To begin with, we estimate the deviation from 1 for the factor f+ defined in Eq. (12). With 
Me = 0.51099828 MeV, w, = 0.114426901 eV, and we ranges from 2 to 7 MeV, the corresponding 


deviation is within the scope +(0.94 — 1.75) %o. Therefore, it is accurate enough to approximate 


f+ as 1. Then, according to the values listed in Table 1, the relative errors of me and w are 


three orders of magnitude lower than that of w2, which means the leading contribution for the 


uncertainty of BCS is the first term in Eq. (11), that is to say we have the relation 


ó&1 = 1 w2 
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with fairly high accuracy (the additional uncertainty is much less than 1078). 

Equation (19) indicates that the feature of electron (positron) beam (denoted by £1 and 
0£1) is totally determined by that of backscattering photon (denoted by w2 and dw2), and vice 
versa. With use of BEMS, w2 is determined by the position of Compton edge while dw2 by the 
slope of the edge!. Herein, the existent of edge slope is just due to óc, the energy spread of 
accelerator. It is noticeable en passant that c; and óc, are constants during measurement. For 
the certain beam energy (fixed £1), the energy spread (6&1) is solely determined by the structure 
of accelerator itself and therefore must be fixed in common sense |15]. However, during the 
data taking of J/V sample performed in year 2012, a peculiar phenomenon is found. As shown 
in Fig. 2, it is noticed that the cross sections vary with the decrease of beam current, which 
means some variations of beam status during energy measurement ?. Such variations imply the 
variant energy, or energy spread, or both of them. The effects, as being elucidated in the next 
section, can lead to energy shift in measurement of beam energy. This is a crux matter for the 


uncertainty analysis of BEMS. 


5 Effects of variations of energy spread and energy drift on mea- 


sured energy 


'This section is devoted to the investigation of effects of variations of energy spread and drift on 
measured energy. Two cases are considered: 1. Assuming there is no energy drift, the energy 
shift is only due to the variation of energy spread; 2. Assuming energy spread is fixed, the 
energy shift is only due to the energy drift. Monte Carlo simulation approach is adopted for the 
following study. 


5.1 Simulation of Compton edge 


The backscattering photons from head-on collision with electron (positron) beam will form a 
sharp edge in a detective spectrum. The pure sharp edge at certain energy (denoted by w) is 


approximated by the normalized function [16] 
h(x) = [p3 + po(x — po)]O(po — 2) . (20) 


Here, the product po(x — po) is small compared to p3; for po = 0 and p3 = 1, h(x) reduces to the 


normal step function. The function h(x) is then folded with a Gaussian of standard deviation 


l'The details of measurement of ws and dw. are delineated in the next section. 
?For BEPCII, in order to control the fluctuation of beam current, a feedback system is added which may affect 


the stability of energy spread. This means the energy spread may change for the fixed beam energy. 
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Figure 2: The relation between beam current and the “relative” cross section with data taken 
at J/w resonance. The horizontal scale is denoted by “Ie x I,/A? ” which means the product of 
electron current (Ie) and positron current (Ip), whose unit is the square of Ampere. The ratio 
of the number of inclusive hadronic events to that of the two-gamma events (both from online 
database), which is proportional to the observed cross section, is denoted as “relative” cross 


section. 
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The resulting function for variable position and height of the edge is given by 
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Anyway, due to the existence of background, a linear function p4(x — po) + ps is added to f(x) 
to describe the shape of background. Therefore, the final synthetic function has the form: 
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The parameters in Eq. (23) are: po - edge position; pı - edge width; po - slope left; p3 - edge 
amplitude; p4 - slope right; ps - background. Parameter po gives the information about the 
average electron beam energy during the data acquisition period, while pı is mostly coupled 
with the electron beam energy spread. 


In the simulation, the following form is adopted 


1 L-W qa ag r-w\? 
f(z) = Zola + a (e = w) ente | c3 slab). (24) 


where w is the position of Compton edge that is used to determine the beam energy; ø is the 


edge width that is related with the beam energy spread (cs); q is unity parameter that is used 
to determine the unity of x (q = 1, the unity of x is MeV while q = 1000, the unity of x is keV). 


Table 2: Input parameters for Compton edge simulation. 


Parameter Value 


Me 0.51099828 MeV 

wo 0.114426901 eV 

Eem 3096.916 MeV 

A 1 MeV 

€ 1548.418 MeV 

Os 0.707107 MeV 
4190.521 keV 

c 3.8272 keV 
1000 

ai —0.1 

ao 0.3 

by —0.2 

bg 0.1 


The relation between w and e is as follows ? 


and the relation between c, and ø is as follows 4 


tJ 
g—2-.—-.0,5. 
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In addition, it should be noted that € = Eem/2 and e, = A/V2, where Eem is the center-of-mass 
(C.M.) energy, and A the spread of Eem. For the energy at J/1» resonance, the input parameters 


for Compton edge simulation are tabulated in Table 2. 


3The relation is just Eqs. (6) and (8), with the correspondence w = w2, € = €1, Wo = w1 and w = w2. 
“The relation is just Eq. (19), with the correspondence € = £1, 75 = 6€1, w = w2 and o = dw. 


5.2 Relation between observed cross section and energy spread 
The cross section of the process e* e^ — J/y — f (where f denotes some final state) is described 
by the Breit-Wigner formula 


127 - Del's 
s— M?)2+T?M?2’ 


oBw(s) = (25) 


where ys is the C.M. energy (s = Eom), Fe and T; are the widths of J/i decaying into 
ete” and f, IT, and M are the total width and mass of J/w. Taking the initial state radiative 


correction into consideration, the cross section becomes [18] 


Orc (W) = [are s) T- Te — pp Gu — z)), (26) 
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wheretm = 1 — s'/s, Vs! is the experimentally required minimum invariant mass of the final 
state f after losing energy due to multi-photon emission; F (s,s) has been calculated in many 
references [18, 19, 20] and II(s(1— zx)) is the vacuum polarization factor. The e*e~ colliders have 
finite energy spread. The energy spread function G(4/s, Vs! ) is usually a Gaussian distribution : 
a 


G(v's, Vs!) = 


27 
vV TA. aq 
where A describes the C.M. energy spread of the accelerator, \/s and Vs’ are the nominal 
and actual C. M. energy respectively. So the experimentally measured resonance cross section 
(observed cross section) is the radiatively corrected Breit-Wigner cross section folded with the 


energy spread function: 


eap(v/s N= [ «(ne Vs! /s)dv/ s! . (28) 


where o... is defined by Eq. (26). 

Numerical calculation indicates that the radiative correction reduces the maximum cross 
section of J/w by 52%; the energy spread further lowers down the cross section by an order of 
magnitude depending on the value of the energy spread. Both the radiative correction and the 
energy spread shifts the maximum height of resonance peak to above the resonance nominal 
mass. In actual experiments, data are naturally taken at the energy which yields the maximum 
inclusive hadron cross section. When the energy spread change, both the maximum cross section 


and the position of energy for the maximum cross section change correspondingly. 


5.3 Relation of cross section with energy spread and energy drift 


The minimization subroutine program DMINFC from CERNLIB [21] is used to find the position 
of energy (Emax) for the maximum cross section and the corresponding maximum cross section 


(Om) itself corresponding to distinctive energy spread (A). The energy shift (Emn) is defined as 
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Figure 3: The relation between the energy spread A and the maximum observed cross section 


Om (a), and the energy shift 0E and the observed cross section ø (b). 
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(b) Fit of Ra 


Figure 4: The relation between beam current and the “relative” cross section, Rọ (denoted by 


diamond in (a)) and “relative” energy spread Ra (denoted by solid circle in (a)). The horizontal 


scale is denoted by “Ie x I,/A? ” which means the product of electron current (Te) and positron 


current (Ip), whose unit is the square of Ampere. 
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(b) Fit of dg 


Figure 5: The relation between beam current and the “relative” cross section, R, (denoted by 


diamond in (a)) and energy drift óg (denoted by solid circle in (a)). 


'The horizontal scale is 


denoted by “Ie x 15/4? ” which means the product of electron current (Te) and positron current 


(Ip), whose unit is the square of Ampere. 
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the difference between the maximum energy (Emax) and the nominal energy (Enom = 3096.916 
MeV), that is 


Emn = Umaz — Enom , 


and the fit curve for 0 E,4,, against A is 
foe, (2) [ke V] = 1.1117 - z[MeV] + 98.567 , (29) 
and the fit curve for om against A is (refer to Fig. 3(a)) 
fon (z)[pb-!] = 3.3885 /x°-8232? — 0.24237 ,  [r: MeV] . (30) 


In Fig. 4(a), the ordinate is the ratio of the number of inclusive hadronic events to that 
of the two-gamma events, which is proportional to the observed cross section and denoted as 
“relative” cross section (R,) in this monograph; in Fig. 4(b), the ordinate is the ratio of two 


energy spread, denoted as “relative” energy spread (RA). The fit curve for Ra against Ie x Ip is 
fra (x) = 0.96451 + 4.5012 - x — 46.612 - z? + 129.44 . 2° . (31) 


Figure 3(b) show the relation between observed cross section e and energy drift 6E, which 
is defined as follows 
ôE = E — Emas , 


where Emar = 3096.9981 MeV corresponding to the energy spread 1 MeV. The fit curve is 
fol(x)[pb7t] = 3.1391 - RM [r: MeV] . (32) 


The ordinate of Fig. 5(a) is the same as that of Fig. 4(a); the ordinate of Fig. 5(b) is the 
energy drift, denoted as (ôg). The fit curve for dz against I. x Ip is 


fa, (x) = —0.38548 + 23.955 - x — 216.47 - £? + 561.15 à? . (33) 


As a matter of fact, in order to obtain the relation between Rọ and Ra (dg), the special 
normalization is adopted. The nitty-gritty is elaborated in Appendix B. Here, fr, (a) in Eq. (31) 
and f5,, (x) Eq. (33) reflect the possible variation of energy spread or energy drift within one beam 
injection, and can be treated as probability function for variation of Ra and dz. Therefore, the 
variable z should take the value between 0 and 1. In the following simulation, the acceptance- 


rejection technique [22, 23] is adopted for distribution sampling. 
5.4 Fitted energy for different cases 


5.4.1 Effect due to variant energy spread 


The simulation for Compton edge is performed for two cases: 1) for the fix energy spread, the 
sampling of the edge is according to distribution formulated in Eq. (24); 2) for the variate energy 
spread, the sampling of the edge is also according to distribution formulated in Eq. (24) but 


with variate c, that is the fix value of ø is replaced by variant c (z), i.e. 
a(x) — o: fra (2), 
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Figure 6: 'T'he comparison of simulations of Compton edge for different cases: fix and variant 
energy spread (a); fix and variant energy drift (b). The fit results for different cases: fix energy 
(c) and variant energy spread (d); fix energy (e), with energy drift (f). Herein, Eem is distinctive 


for (c) and (d) as indicted in Table 4. In all plot circle denote the fix cases while box the variant 
cases; the curve is the best fit result. There are 0.5 Million counts in each sample. 
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where fr, (x) is the distribution in Eq. (31), and x is a random number of between 0 and 1. 
The simulated distributions for two cases are shown in Fig. 6(a). 

The fitted results based on the simulated distributions for two cases are given in Table. 3, 
and the fit curves for two cases is displayed in Fig. 6(c) and (d). The difference for fitted beam 
energy is about 0.113 MeV, or 0.226 MeV for C.M. energy. 


Table 3: The fit results of beam energy for different cases of energy spread : fix energy spread 
and variant energy spread. 


Parameter Fix a Variant c 

x? 83.64406 77.50543 

€ 1548.300 + 0.082 MeV 1548.413+ 0.115 MeV 
A 1.1215 + 0.0480 MeV 1.0924 + 0.0573 MeV 
ài —0.12464 + 0.01937 —0.10374 + 0.03012 

ag 5.8016 + 0.0281 5.7938 + 0.0285 

by —2.0528 + 0.0807 —2.0300 + 0.0819 

bg 1.9456 + 0.0134 1.9471 + 0.0133 


Table 4: The fit results of beam energy for different cases of energy drift : fix energy and with 
energy drift. The input Eem = 3096.9981 GeV instead of Eem = 3096.916 GeV. 


Parameter Fix w Variant w 

y 86.93219 76.01074 

€ 1548.339 + 0.078 MeV 1548.572 + 0.100 MeV 
A 1.1211 + 0.0469 MeV 1.0631 + 0.0517 MeV 
ay —0.12588 + 0.01884 —0.11371 + 0.02701 

ag 5.7998 + 0.0280 5.7678 + 0.0283 

by —2.0572 + 0.0804 —2.0276 + 0.0822 

b» 1.9461 + 0.0134 1.9384 + 0.0132 


5.4.2 Effect due to energy drift 


The simulation for Compton edge is performed for two cases: 1) for the fix energy, the sampling 
of the edge is according to distribution formulated in Eq. (24); 2) for the variant energy, the 
sampling of the edge is also according to distribution formulated in Eq. (24) but with variant 


w, that is the fix value of w is replaced by variate w(x), i.e. 


w(x) =W: fy (2) , 


where f5,, (x) is the distribution in Eq. (33), and x is a random number of between 0 and 1. The 
simulated distributions for two cases are shown in Fig. 6(b). 

The fitted results based on the simulated distributions for two cases are given in Table. 4, 
and the fit curves for two cases is displayed in Fig. 6(e) and (f). The difference for fitted beam 
energy is about 0.233 MeV, or 0.466 MeV for C.M. energy. 
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6 Summary 


In this monograph, the energy relation between Compton backscattering photon and high en- 
ergy electron is derived analytically, based on which the formula for uncertainty estimation of 
measured energy is obtained. The leading contribution of uncertainty is figured out by utilizing 
the present experimental information. 

By virtue of the experimentally available information, the special phenomenon between beam 
current and cross section is explored in detail by simulation approach. The quantitative results 
indicate that for C.M. energy, the maximum energy shift is up to 0.226 MeV due to variant 
energy spread and 0.466 MeV due to variant energy drift. These effects of energy shift studied 
herein disclose so to speak a significantly possible source of systematic uncertainty for BEMS, 


which in turn has the far-reaching meaning for the further analysis of physics error. 
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Appendix A 


The nature unity with c = 1 is adopted hereafter. The velocity of electron is denoted as 3(= v/c) 
and its corresponding energy is often expressed as € = my with y = 1/4/1 — 8?. According to 
Refs. [24, 25], the general relation between w; and wz for Compton scattering process is 


T «1(1 — f cos $1) (A.1) 


tapot T ceased dal) 
ym 


where $4 is the angle between incident photon and electron while $5 the angle between backscat- 


tering photon and electron. For head-on collison, $1 = 180° and $9 = 0°, so Eq. (A.1) becomes 


1+8 
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Note two relations 
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and also the relation m?4? = £2, Eq. (A.2) can be rewritten as 
er 8)? 
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and this is just Eq. (5), from which it yields 
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Then the partial derivative with respect to w1, w2, and m can be expressed by 
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Based on the law of error propagation, the 9&1 is obtained as follows 


(861)? = (ze 7 IE (= IE (m) (A.7) 


Appendix B 


In this appendix, c, 6 and x represent respectively the relative cross section (Ro), energy spread 
(RA), and energy drift (0E). Firstly, 


gj. = Dg +r, 
(A.8) 
quu = z + pi. 
Some algebra yields 
Do BUE. 
r po+ ka 
with definitions 
Pr, kl. a; = pı - 0," eq. 
O1 ay 
Then it is easy to acquire 
— — (A.9) 


(p2 --ai)r — ay 

With the above relation, if c4 and 6, (equivalently a4) are chosen, à» can be calculated from o2». 

In our study, c1 is chosen as the maximum cross section, which guarantees k is always greater 

than 1. As to 6), 0.8 MeV, 1.0 MeV, and 1.2 MeV, are used to calculate the o», the shape of 
curve are exactly the same. Without loss of generality and actuality, 6; is set to be 1.0 MeV. 
Secondly, we consider the relation between o and z. 

cy = pne, 


0 = pe P, iy) 
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Some algebra yields 
K 
— = n(n + 221) , 
pa 


with definitions 


01 
K=In—, mg-z—m. 
02 


According to root formula for quadratic equation and notice 7 > 0, then 


£2 = 4/22 + K/p2 . (A.11) 


Let xı correspond to c4 and e, to the maximum cross section, that is xı = 0 and 64 = Omaz, 


then we have 
1 Omar 
-In —— . 


p2 02 


(A.12) 
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